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A combination of parabolized stability equations and secondary instability theory has been 
applied to a low-speed swept airfoil model with a chord Reynolds number of 7.15 million, 
with the goals of (i) evaluating this methodology in the context of transition prediction for a 
known configuration for which roughness based crossflow transition control has been 
demonstrated under flight conditions and (ii) of analyzing the mechanism of transition delay 
via the introduction of discrete roughness elements (DRE). Roughness based transition 
control involves controlled seeding of suitable, suhdominant crossflow modes, so as to 
weaken the growth of naturally occurring, linearly more unstable crossflow modes. 
Therefore, a synthesis of receptivity, linear and nonlinear growth of stationary crossflow 
disturbances, and the ensuing development of high frequency secondary instabilities is 
desirable to understand the experimentally observed transition behavior. With further 
validation, such higher fidelity prediction methodology could be utilized to assess the 
potential for crossflow transition control at even higher Reynolds numbers, where 
experimental data is currently unavailable. 


I. Introduction 

S ubstantial reductions in the fuel burn of future aircraft have been targeted both in the U.S. [1] and in Europe 
[2], in order to lower the cost of air travel and alleviate the impact of aviation on the environment. A 
substantial fraction of the targeted reductions in aircraft fuel burn will accrue from the reductions in vehicle drag. 
For today’s commercial transport aircraft, up to one half of the total drag corresponds to skin friction drag. Because 
laminar skin friction is much lower in comparison with the turbulent skin friction, flow control via delayed 
boundary-layer transition over major aerodynamic surfaces holds the potential to provide the desired reductions in 
the overall drag. 

NASA’s Subsonic Fixed Wing Project is investigating several options to help develop robust technology 
concepts for laminar flow control and to accelerate their insertion into operational fleet [1]. One of those concepts 
involves the delay of crossflow-induced transition over a swept aerodynamic surface using control input in the form 
of appropriately positioned Discrete Roughness Elements (DREs) near the leading edge [3]. A bulk of the existing 
experimental and computational studies of the DRE concept have been carried out for either incompressible or 
supersonic configurations with modest wing-chord Reynolds numbers of up to approximately 8 million [4, 5], and 
with pressure distributions that may or may not be optimal for wing designs for subsonic transport aircraft flying at 
Mach numbers between 0.75 to 0.90. Consequently, additional effort will be necessary to assess the potential 
capability of the DRE concept to control swept wing transition at transonic Mach numbers and substantially higher 
chord Reynolds numbers than previous applications. 

Given the sensitivity of transition to wind tunnel disturbances and the cost of flight experiments, computational 
tools will need to play an important role in the assessment of DRE for subsonic transport vehicles, as a means of risk 
reduction as well as to help optimize the DRE design. The DRE concept seeks to delay transition via controlled 
seeding of subdominant crossflow modes that cannot lead to transition on their own, but can keep the naturally 
dominant instability modes at bay via nonlinear modification of the basic state and/or mode competition. This 
intrinsically nonlinear control mechanism calls for a higher fidelity prediction approach involving most, if not all, of 
the relevant stages of the transition process, namely, (1) excitation of primary crossflow instability, (2) linear and 
nonlinear development of the primary crossflow instability, and (3) excitation, linear growth and nonlinear 
breakdown of high-frequency secondary instability [6] . A number of individual ingredients of this approach have 
been available for some time and partial integration of those ingredients has been demonstrated in prior works 
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[6-11]. The important missing parts are the receptivity and the nonlinear breakdown of secondary instability, the 
latter of which directly leads to the onset of transition as manifested by a sudden rise in the wall shear-stress. 

The objective of the present work is to utilize a combination of parabolized stability equations and secondary 
instability theory, to (1) assess the effectiveness of DRE control and (2) compute the nonlinear development and 
breakdown of secondary instability up to the onset of transition. The capability to carry out the latter task implies 
that the high fidelity, fully integrated transition prediction tool is one step closer to reality. 

The receptivity of the stationary crossflow vortex and that of the secondary instability will not be discussed in 
this paper. A number of investigations in existing literature have addressed the former problem ( e.g . [9, 11]), but the 
latter problem does not appear to have been tackled so far. Thus, parametric studies encompassing a range of initial 
amplitudes of the stationary crossflow vortex as well as the secondary instability modes are used to gauge the 
resulting variations in transition characetristics. The rest of the paper is laid out as follows. Section II briefly 
describes the flow configuration under examination and the computational codes used during the present effort. In 
Section III, the results of stationary crossflow computations are presented and the effect of control mode on target 
mode evolution is analyzed. In section IV, results of linear secondary instability computations are presented and the 
effects of control on linear growth characteristics of the secondary modes are discussed. In section V, PSE 
computations of the nonlinear development of secondary instability modes, the ensuing breakdown process and 
transition criteria accounting for the secondary instability characteristics are discussed. Finally, a summary of the 
paper is given in Section VI. 


II. Flow Configuration and Analysis Codes 

The flow configuration examined herein corresponds to the SWIFT flight test article from [4]. A 30-degree 
swept airfoil was mounted vertically beneath the wing of a Cessna 0-2A Skymaster aircraft and flown at chord 
Reynolds numbers of 7 to 8 million, in order to investigate the feasibility of DRE technology for SensorCraft 
applications. The mean flow around the installed configuration was computed via CFD [5]. Span wise variations in 
the computed pressure distribution over the mid-span segment of the suction surface were relatively small, so that 
the region in between the rows of static pressure ports at 13 inch and 29 inch span locations outboard of the 
mounting plane approximated (locally) infinite span conditions. Transition measurements, both without and with 
the use of DRE, were reported in [4] . 

All of the results presented in this paper correspond to the inboard measurement location at Test Point 27 [4]. 
For this case, the Reynolds number Re c based on the free-stream velocity and the streamwise chord length is equal 
to 7.15xl0 6 and the angle of attack over the test article is -4.69 degrees. The mean flow at these conditions 
computed with a boundary layer solver [12] using infinite span approximation, in conjunction with the sectional 
surface pressure distribution derived from the CFD solution [5]. Linear and nonlinear development of the instability 
modes is computed using parabolized stability equations (PSE) as implemented in the Langley Stability and 
Transition Analysis Codes (LASTRAC) [13]. Linear secondary instability predictions are made using the 
methodology outlined in [10], whereby the growth of high frequency instability modes supported by the finite 
amplitude stationary crossflow vortex is analyzed in a manner similar to classical linear stability analysis of swept 
airfoil boundary layers. The main difference between the primary and secondary stability analyses is related to the 
more complex basic state for the secondary instability analysis. Specifically, this basic state consists of the mean 
boundary layer flow together with the finite amplitude primary crossflow mode and, therefore, varies in both surface 
normal and spanwise direction directions. Consequently, the instability characteristics of the secondary modes must 
be analyzed using a planar, partial differential equation based eigenvalue problem, rather than as an ordinary 
differential equation based eigenvalue problem for the classical analysis. Nonlinear development and breakdown of 
the secondary instability are again computed using LASTRAC. The selection of grid parameters and other aspects of 
the numerical solution are based on extensive experience with such calculations for similar classes of flows [6, 7, 9, 
10], and spot checks were made to ensure that any variations with respect to modifications to those choices were 
sufficiently small. 


III. Stationary Crossflow Vortices and Effect of Control 

Linear amplification characteristics of stationary crossflow instability at the flow conditions of interest have been 
described in [5]. Specifically, the crossflow mode with a spanwise wavelength of X = 4.5 mm corresponds to the 
maximum logarithmic amplification ratio, N - 12, based on linear stability theory. Nearby waves of wavelengths 
between 4 and 5 mm have comparably large amplification ratios. The first harmonic of the 4.5 mm mode (k = 2.25 


2 



mm) has a considerably lower linear growth potential (N ~ 5.6) and, therefore, was used in the experiment [4] to 
control the growth of the more amplified instability modes via DRE. In what follows, some of the computations are 
carried out for the 5 mm vortex instead of the 4.5 mm mode. Since the linear amplification ratios for both modes are 
close to each other, no major differences in computational results are expected. For control action computations, 
however, the 2.25 mm and the 4.5 mm modes are consistently used. 

Flight observations under somewhat different conditions (Re c = 8.1xl0 6 and AO A = -4 deg) [4] had shown that, 
with a painted leading edge having an average roughness height of 1.0 pm rms and 3.8 pm average peak-to-peak, the 
natural onset of transition at the span wise section of interest was near x/c ~ 0.3, whereas transition was delayed 
beyond x/c ~ 0.8 when the leading edge surface was unpainted and, therefore, was considerably smoother with a 
r.m.s. roughness height of 0.3 pm rms and 2.2 pm average peak-to-peak. Even with the painted leading edge, 
transition onset was pushed back to x/c ~ 0.6 when a 12 pm height roughness array with 2.25 mm spacing was 
introduced just downstream of the airfoil leading edge. 

It is reasonable to assume that both the painted and polished leading edges of the test article had a broadband 
surface roughness distribution, with a wavenumber spectrum that overlapped with the previously discussed range of 
most unstable stationary crossflow modes. A match in length scales would lead to a direct excitation of a 
correspondingly broad spectrum of unstable modes via a linear receptivity mechanism [14]. Of course, following 
sufficiently strong and differential amplification over a finite downstream distance, the disturbance spectrum would 
eventually tend be dominated by the modes in the vicinity of the linearly most amplified wavelength. As a simple, 
conceptual model for the effects of background surface roughness on transition via the dominant range of 
wavelengths, the streamwise evolution of a stationary crossflow mode is first considered with X = 5 mm and a range 
of initial amplitudes (Figure la). As previously mentioned, the span wise wavelength of the selected mode (X = 5 
mm) is very close to the wavelength of the linearly most unstable mode corresponding to X = 4.5 mm. Analogous 
evolution of crossflow modes with shorter and longer wavelengths, respectively, is shown in Figures, lb and lc, 
respectively. The limited range of chordwise locations included in each plot is due to the beginning of a strongly 
adverse gradient slightly downstream of x/c = 0.70, which leads to flow separation within the framework of 
boundary layer equations. It is interesting to observe that, regardless of the disturbance wavelength, crossflow 
disturbances with sufficiently large initial amplitude exhibit a peak in the fundamental amplitude, followed by a 
significant decrease in fundamental amplitude in the immediately downstream region. In other words, no extended 
regions of quasi- saturation observed for the lower Reynolds number, NLF-0415(b) configuration [3, 7] are noted in 
the present case. 

Figure 2 illustrates a typical sequence of velocity contours during the nonlinear evolution of an unstable 
crossflow mode. The initially undulating contours give way to strongly overturned contours that are symptomatic of 
strongly inflexional velocity profiles in the outer part of the boundary layer and, hence, the onset of high-frequency 
secondary instabilities that typically, but not always, harbinger the onset of laminar breakdown to turbulence. 

Now, the control scenario is considered. The X = 2.25 mm mode with a larger initial amplitude is introduced as 
the control input for its more unstable subharmonic with X = 4.5 mm. The term “amplitude” is defined here as the 
maximum perturbation velocity in the normal chord direction (i.e. parallel to the leading edge) as a fraction of the 
freestream velocity. The combined evolution of the control and target modes with the initial amplitude of the 4.5 

mm mode held fixed at lO -4 (based on the preceding discussion for the 5 mm mode) is shown in Figure 3(a). As the 
amplitude of the control input is increased from 0 to 0.01, the amplitude levels of the target mode within the 
downstream region become substantially smaller. Of course, the peak amplitudes of the control mode also increase 
with the amplitude of the initial control input, which introduces the possibility of the control mode sustaining high- 
frequency secondary instabilities and, hence, precipitating an earlier onset of transition. This possibility will be 
discussed in the subsequent sections of this paper. 

As shown in Figure 4(a), at least part of the stabilization may perhaps be attributed to the reduced crossflow 
Reynolds numbers in the presence of the control mode with nonlinear amplitudes [15]. However, the same control 
input appears to result in larger values of R^ farther downstream; therefore, the mechanism underlying the 
stabilization of the target mode cannot be attributed to changes in R^ alone. Figure 4(b) compares the linear N- 
factors of the target mode as predicted via instability analysis of the spanwise averaged mean flow without and with 
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the control input. The figure shows that the modified basic state leads to an N-factor for the target mode that is 
smaller in comparison with the unperturbed flow by approximately 3.2. 

IV. Linear Secondary Instability and Effect of Control 

For stationary crossflow vortex instability, the dominant velocity component is that along the vortex axis and the 
velocity components in the cross-plane normal to the vortex axis are relatively small in comparison. In some of the 
previous investigations ( e.g . [16]), the flow within the cross-plane has, nonetheless, been found to have a large 
influence on the growth rates of high frequency secondary instability supported by the finite amplitude primary 
crossflow modes. To assess the importance of the cross-plane flow, the secondary instability characteristics were 
computed both with and without the cross-plane velocity components associated with the distorted mean flow. The 
secondary instability growth rates at x/c ~ 0.45 are shown in Figure 5(a) for a primary crossflow vortex with X = 5.0 
mm and an initial amplitude of lO -4 . Clearly, those components appear to have a rather small effect in this 
particular case. It is found that there are multiple families of unstable secondary modes, but only the three most 
dominant ones are shown in Figure 5(a). Representative mode shapes for these dominant modes are shown in 
Figure 5(b). The included labels for each family (F, Z, and Y/Z) indicate whether the wall-normal and/or the 
spanwise shear associated with the distorted basic state (Figure 2) dominates the energy production mechanism 
associated with a given family of secondary instability modes [7]. The FZ-mode looks like a F-mode in the sense 
that it rides on top of the stationary crossflow vortex with a mode shape similar to that of a F-mode (see bottom of 
Figure 2(b)); however, it is categorized as a mixed Y/Z - mode because its energy balance reveals that the contribution 
to its growth rate from the production term associated with the spanwise shear is comparable to the contribution due 
to wall-normal shear. 

The spatial amplification characteristics of various fixed frequency secondary instabilities belonging to each of 
the three dominant families of secondary instability modes are shown in Figure 6. It is observed that the secondary 
modes grow extremely rapidly in comparison with the finite amplitude primary crossflow modes, achieving 
(secondary) N factors in excess of 9 across a short distance of around 6 percent chord. Of course, the receptivity 
characteristics of each of these modes will determine whether or not that mode can get excited with an adequately 
large initial amplitude to have a visible influence on the transition process. For example, previous experimental 
studies [17, 18] indicate a preferential excitation of the Z-modes under certain conditions. Given this scenario, it 
seems prudent to track the N-factor envelopes for each mode family separately rather than tracking a combined 
envelope for all secondary instability modes. 

The effect of control on primary stationary crossflow vortices has been discussed in the previous section. The 
associated evolution of secondary instability is studied here. With the initial amplitude of the 4.5 mm primary 

stationary crossflow vortex mode fixed atlO -4 , the effect of increasing the amplitude of the 2.25 mm control mode 
is shown in Figure 7. With zero control amplitude, it is seen that both the dominant F-mode and Z-mode undergo 
strong modal amplifications and that the difference in estimated transition onset location would be Ax/c = 0(0.05) 
depending on whether the F or the Z modes are chosen for transition correlation. 

As the control amplitude increases, modal amplifications of both the dominant F- and Z-modes are reduced. As 
shown in Figure 7, the presence of the control mode apparently affects the Z-mode of secondary instability much 
more than the F-mode. This indicates that, if indeed transition is more often caused by the growth of the Z-mode as 
some experiments seem to suggest [17, 18], then DRE is a very effective means to achieve transition delay in the 
present case. 

Note that the mixed Y/Z mode has been conspicuously missing from the above discussion of the effect of control 
input. Without control, the Y/Z mode is subdominant to the F-mode. With increasing control amplitudes, this mode 
is also suppressed. An examination of its energy balance shows that the contribution from the part of the production 
term associated with the spanwise shear becomes negligible or even negative, i.e., the mixed mode now transforms 
into a predominantly F-mode and is always subdominant to the most unstable F-mode shown in Figure 7. For all 
practical purposes, no special attention needs to be paid to this mixed Y/Z mode, except to note that the morphing 
from a mixed Y/Z - mode to a pure F-mode is further evidence that, at the present set of flow conditions, the DRE is 
more effective in suppressing the instability mechanisms associated with the spanwise shear. 
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While introducing the subdominant crossflow mode at relatively larger initial amplitudes may lead to substantial 
reduction in the growth of the otherwise dominant crossflow mode and, hence, in the growth potential of the 
secondary instability modes associated with the latter mode, the possibility of overcontrol (i.e., risk of premature 
transition due to the finite amplitude control mode) cannot be eliminated. That is, in this case, the 2.25 mm control 
mode could become susceptible to secondary instability and may itself break down to cause an earlier onset of 
transition. Secondary instability analysis for the 2.25 mm control mode is carried out to assess the likelihood of 
such a scenario. The peak N-factors for Z-mode of secondary instability are less than 1 for all three control inputs. 
Of the three initial amplitudes of 0.002, 0.005 and 0.01, only the last gives rise to a maximum F-mode secondary 
instability N-factor of approximately 11 at a chordwise location close to x/c = 0.25. For the other two amplitudes of 
control input, the F-mode N-factors remain less than or equal to approximately 8. Therefore, even if the F-mode can 
also lead to transition, the possibility of premature transition due to over control is remote as long as the control 
amplitudes are kept under a certain threshold, which in this case appears to be around 0.005 or slightly larger. 

V. Nonlinear Development and Breakdown of Secondary Instability 

As described in section IV, the linear predictions of secondary instability may be used for transition onset 
correlation using a secondary N-factor criterion. However, it is also instructive to study the nonlinear evolution of 
the secondary mode perturbations in order to study the breakdown mechanisms responsible for the onset of 
transition. Wassermann and Kloker [19] carried out direct numerical simulations (DNS) of secondary instability 
well into the breakdown region and revealed some interesting features of the breakdown process for a specific wind 
tunnel configuration. Although DNS is a powerful tool, it is also costly. The present section investigates the 
application of nonlinear parabolized stability equations (PSE) to capture the initial portion of the breakdown process 
at a computational cost that is at least an order of magnitude lower in comparison with the DNS. Illustrative PSE 
computations of the breakdown of crossflow vortices with X = 5 mm via the secondary instability mechanism are 
presented below. 

The secondary instability eigenfunction is computed at a short distance downstream of the neutral point 
predicted by the linear, 2-D eigenvalue analysis, and is fed into the nonlinear PSE computation as initial conditions 
at that location. Figure 8 shows the comparison between the secondary instability growth obtained via PSE and 
linear theory, respectively. The left half of Figure 8 shows the evolution of a Y-mode of secondary instability at the 
frequency of f = 36 kHz, whereas the right half shows similar comparison for a Z-mode at f = 32 kHz. The red and 
green lines represent, respectively, the amplitudes of stationary crossflow vortices and secondary instability at 
different spanwise wave numbers obtained by PSE computations. The black dots represent the linear growth 
predictions based on the secondary instability theory. It is seen that, when the secondary instability amplitude is 
small, PSE results and the quasi-parallel linear theory agree very well for the Z-mode. For the F-mode, there is some 
initial discrepancy between the linear theory and PSE results. The reason for this is currently unclear, but the PSE 
solution is known for having initial transient behavior when the inflow condition deviates from the exact non- 
parallel eigensolution. In any case, after the transient dies out, the PSE solution again agrees very well with the 
secondary instability theory. 

The transition location can be correlated to either the amplitude of the nonlinearly developing secondary 
instability or the rise in wall-shear stress that characterizes the onset of transition. As discussed before, breakdown 
of crossflow vortices via the Z-mode of secondary instability has been observed in limited measurements and 
computations reported in the literature [16, 17, 18]. While this does not necessarily rule out the potential relevance 
of the F-mode of secondary instability, it appears reasonable to examine the nonlinear evolution of the Z-mode of 
secondary instability first. Thus, a parametric study of the transition onset location as a function of the initial 
amplitude of secondary instability is carried out. Figure 9 shows the normalized wall shear stress (red lines), 
together with the RMS amplitude (green lines) of the secondary instability for several initial amplitudes of the Z- 
mode. The left most shear and amplitude curves correspond to an initial chordwise velocity perturbation of 0.01, and 
successively adjacent curves to the right are associated with lower initial amplitudes in 10-fold decrements. In all 
cases, as the secondary instability amplitudes become sufficiently nonlinear, the wall-shear curve starts to rise and 
becomes nearly vertical at a short distance downstream. This rapid rise in wall shear is symptomatic of the onset of 
transition. If the transition onset location is to be based on the location where the wall-shear stress rises sharply, then 
the uncertainty in the predicted transition onset location (indicated by the dashed vertical in Figure 9) is very small 
provided that the initial amplitude of the secondary instability is known. 
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The transition onset location could also be defined in terms of the location where the root mean square (RMS) 
amplitude of the secondary instability exceeds a certain threshold value. However, this transition criterion may be 
somewhat less robust due to the fact that the RMS amplitude does not rise as sharply as the wall shear and, in some 
cases, also undergoes a drop-and-rise in the vicinity of the transition onset (see green lines in Figure 9). 

It can be further noted in Figure 9 that, with each increase in initial amplitude of the most amplified secondary 
instability mode by a factor of 10, the transition location based on, say, the wall- shear criterion moves upstream by 
approximately equal intervals. The red line in Figure 10 which represents the transition location as a function of the 
logarithm of the initial amplitude of the secondary instability is almost a straight line, at least in this particular case. 
This appears to indicate that the nonlinear effect due to the amplification of higher harmonics is limited to very short 
streamwise domain just upstream of the transition onset. The implication is that linear secondary instability N-factor 
methods, together with known initial amplitude of the secondary perturbation, may be sufficient to predict the 
transition location with reasonable accuracy. Ref. [20] describes linear amplitude criteria based on primary 
instabilities and their application for transition prediction. 

Suppose that the transition onset is correlated with a constant linear amplitude of A tr . Then, by linear theory, 
the relation between the initial amplitude and the N-factor at transition is, 

Ar = A e N,r (1) 

where A. and N tr represent the initial amplitude and the transition N-factor, respectively. Because A tr is assumed to 
be constant, it follows that 

^ tr = _i (2) 

Alog(A) ’ ( } 

i.e., the slope of the N tr vs. log( A. ) curve is equal to -1. 

The linear N-factors at the transition locations resulting from the Z-mode at 33 kHz are also plotted (represented 
by the green line) in Figure 10. The curve is approximately a straight line whose slope (represented by the dashed 
line) is very close to that given by Eq (2). 

Figure 11 shows the RMS amplitude distributions associated with the high frequency instability at selected 
chordwise locations. At each station, the amplitude distribution is normalized so as to have a maximum value of 
unity. For small values of x/c , the amplitude distributions mimic the linear eigenfunctions shown in Figure 3 and 
remain so over much of the region between the locations where the secondary instability is initiated and where the 
transiton onset occurs. Very close to the onset of transition, the amplitude distributions begin to be less concentrated 
and start to approach the wall, eventually leading to the sharp rise in spanwise averaged wall shear. 

In the results presented thus far, the visualization of secondary instability were limited to 1- or 2-dimensional 
representations. The descriptions below are intended to show a rudimentary 3-dimenional picture pertaining to the 
nonlinear growth of secondary instabilities. Without secondary instability, the U-iso-surfaces for purely stationary 
crossflow vortices would look like smooth tubes that are approximately aligned with the boundary-layer edge 
velocity. With the introduction of the secondary instability, these iso-surfaces are modified. The top and bottom 
portions of Figure 12 depict the instantaneous iso-surfaces of chordwise velocity, U, for the cases initiated with 7- 
and Z-modes of secondary instability, respectively. The initial amplitudes of the secondary instability modes for 
both cases in Figure 12 are equal to 0.01. The “rib-like” structures previously observed in [19] and [21] are clearly 
visible in both cases. By examining the vorticity of these “ribs”, it can be determined that these are, in fact, co- 
rotating vortices which ride on the top of the stationary crossflow vortex in the case of 7-mode, or on the side of the 
stationary vortex in the case of Z-mode of secondary instability. Furthermore, these “ribs” align themselves at an 
oblique angle relative to the axis of the stationary crossflow vortex (see, also, ref. [19] in this context). Farther 
downstream, smaller structures begin to appear, which should eventually lead to turbulence. 

VI. Summary 

The present effort targets the development, application, and improvement of higher fidelity tools for DRE based 
laminar flow control over subsonic aircraft. The results obtained for the SWIFT flight configuration demonstrate 
that useful insights may be gained into the efficacy and/or optimization of a given DRE application using such tools. 
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Computations of the nonlinear evolution of stationary crossflow modes using the parabolized stability equations 
reveal that, with increasing amplitude of the control mode with X = 2.25 mm, a corresponding reduction in the 
amplitude of the target mode with X = 4.5 mm is achieved. 

The reduction in the target mode amplitude is accompanied by reduced linear growth rates of the secondary 
instabilities. Several secondary instability modes are found, of which three representative ones with the largest 
growth rates are analyzed. It is revealed that the presence of the control modes reduces the growth of the Z-mode 
secondary instability (i.e., modes primarily driven by disturbance energy production associated with the span wise 
shear of the basic state) much more effectively than the F-modes (that derive a majority of their amplification from 
the wall-normal shear of the basic state). 

The possibility of over-control is also assessed, and it is found that, if and when the F-mode is a catalyst for 
transition, then an initial control amplitude exceeding a certain threshold can make the control mode itself 
susceptible to secondary instability breakdown, leading to an earlier-than-intended transition at the flow conditions 
investigated in this paper. However, if the Z-mode alone is responsible for transition, then this threshold amplitude is 
significantly higher in comparison. 

The nonlinear development of the most unstable Z-mode of the secondary instability is traced using the PSE 
method. The PSE results agree well with linear secondary instability when the amplitude is small. At large 
amplitudes, a sharp rise in the wall-shear stress indicates the onset of transition. 

The role of low frequency, traveling crossflow instabilities during the transition process has not been discussed 
in this paper. It is known that, in the presence of strong stationary crossflow vortices, traveling crossflow 
instabilities are modulated [7, 10], with maximum amplitude nearer to the wall than the secondary instabilities. 
Therefore, future studies may include the development of modulated low frequency, traveling crosssflow 
instabilities and their interaction with high frequency secondary instabilities. 
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Figure 1 . Nonlinear evolution of fundamental disturbance amplitude for a stationary crossflow mode at specified 
spanwise wavelength, with the initial disturbance amplitude as a parameter 



Figure 2. Chordwise velocity contours across vortex cross section extending over one spanwise wavelength 
(k = 5 mm, Aj = le-4) 
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Figure 3. Modal amplitude evolution of linearly most unstable stationary crossflow mode in the presence of 
control input at k = 2.25 mm 



(a) Relative change in crossflow Reynolds number 
distribution as a result of control input at various 
levels 



naturally dominant stationary crossflow mode (k 
= 4.5 mm) with and without the mean flow 
correction due to control input with an initial 
amplitude of Ai = 0.01 


Figure 4. Effect of control input (k = 2.25 mm) on instability characteristics of span wise averaged flow 
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(a) Temporal growth rates with (solid lines) and without 
(dashed lines) V 0 and W 0 


F-mode 


Z-mode 



F/Z- mode 

(b) Mode shapes for three dominant families of modes 


Figure 5. High-frequency secondary instabilities at x/c = 0.45 (primary crossflow mode corresponds to span wise 
wavelength of X = 4.5 mm and initial amplitude of Ai = le-4) 





Figure 6. N-factor curves for spatially growing 
disturbances at selected frequencies within each 
family of high-frequency secondary instabilities 
(indicated by different colors as described in the 
legend). 



Figure 7. Envelopes of N-factor curves for dominant 
Y and Z families of secondary modes at 
increasing levels of control amplitude. 
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Normalized wall shear 



X/C x/c 


Figure 8. Amplitude evolution of primary and secondary instability perturbations. Each line is associated with a 
fixed spanwise wavenumber. The red lines represent stationary modes (predominantly associated with 
primary crossflow instability) and the green lines represent secondary instability modes at the fundamental 
frequency (higher frequency harmonics are not shown). The black dots represent results from linear 
secondary instability theory. The y-mode secondary instability is shown on the left, and the z-mode on the 
right. 



Figure 9. Evolution of wall shear and secondary 
instability RMS amplitude for Z-mode 
secondary instability as control amplitude 
varies. The dashed vertical lines indicate the 
transition locations based on the wall-shear 
criterion. 



Figure 10. Transition location and transition N-factor as 
function of logarithm of secondary instability 
initial amplitude is almost a straight line. The 
latter falls approximately on a straight line with 
a slope given by Eq (2). (Note: the expression 
indicating the slope is equivalent to Eq (2), but 
is expressed in logarithm to the base 10 instead 
of natural logarithm to be consistent with the 
scaling of the horizontal axis). 


12 



Figure 11. Time-RMS distributions of chordwise velocity distributions at selected chordwise stations. The white lines 
in the background are the u-contours of stationary crossflow vortex. The y-mode is displayed on the left and 
the z-mode on the right. Secondary instability initial amplitude is le-6. 
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Figure 12. Iso-surfaces of instantaneous chordwise velocity U. Top: F-mode, U =0.91t/ oo . Bottom: Z-mode, 
U =0.87£/ oo 
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